Driven interfaces in random media at finite temperature : 
is there an anomalous zero-velocity phase at small external force ? 
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The motion of driven interfaces in random media at finite temperature T and small external force 
F is usually described by a linear displacement hdt) ~ V(F,T)t at large times, where the velocity 
vanishes according to the creep formula as V(F,T) ~ e - K ( T )/ F f or p fj. In this paper, we 
question this picture on the specific example of the directed polymer in a two dimensional random 
medium. We have recently shown (C. Monthus and T. Garel, arXiv:0802.2502 ) that its dynamics 
for F = can be analyzed in terms of a strong disorder renormalization procedure, where the 
distribution of renormalized barriers flows towards some "infinite disorder fixed point". In the 
£N| ' present paper, we obtain that for small F, this "infinite disorder fixed point" becomes a "strong 

' disorder fixed point" with an exponential distribution of renormalized barriers. The corresponding 

distribution of trapping times then only decays as a power-law P(t) ~ 1/r , where the exponent 
a(F,T) vanishes as a(F,T) ex _F M as F — > 0. Our conclusion is that in the small force region 
a(F, T) < 1, the divergence of the averaged trapping time r = +oo induces strong non-self-averaging 
effects that invalidate the usual creep formula obtained by replacing all trapping times by the 
typical value. We find instead that the motion is only sub-linearly in time ha(t) ~ t a ^ F ' T \ i.e. 
the asymptotic velocity vanishes V = 0. This analysis is confirmed by numerical simulations of a 
directed polymer with a metric constraint driven in a traps landscape. We moreover obtain that the 
I . roughness exponent, which is governed by the equilibrium value ( eq = 2/3 up to some large scale, 

, ^ ' becomes equal to £ = 1 at the largest scales. 

: 

I. INTRODUCTION 

0: 

| A. Dynamical phase diagrams in the presence of quenched disorder 

c : 

Transport phenomena in random media have remained a very active field of research since the discovery of the 
Anderson localization fifty years ago For classical systems also, the presence of quenched disorder can induce 
completely new transport behaviors with respect to pure systems in some regions of parameters. It turns out that 
(T) ■ even the dynamics of a single particle in one-dimensional random media (see the review Q and references therein) 
can already present a very rich phase diagram as a function of the temperature T and the external applied force F. 
As an example, we shown on Fig. 1 the exactly known phase diagram for the biased Sinai model [J 0, d, @ : one 
needs to introduce a dimensionless parameter a(F,T) — TF (where we have chosen to fix the disorder strength at 
some simple value to simplify the notations and emphasize the dependence upon T and F we are interested in ) . For 
F = 0, the motion is logarithmically slow x(t) ~ ±(lnt) 2 . For < a(F,T) < 1, the diffusion is anomalous with a 
sub-linear displacement x ~ £ Q ( F < T ) ■ f or 1 < a < 2, the velocity becomes finite x ~ Vt, but the dispersion remains 
anomalous of order t 1 / ct ( F ' T )_ Finally for a > 2, both the velocity and the diffusion coefficient are finite. So already 
in simple one-dimensional systems, various dynamical properties can undergo phase transitions at various thresholds. 

The case of more complex systems like interfaces or manifolds in random media has attracted a lot of attention in 
relation with many important applications, and various approaches have been developed to elucidate the structure 
of the dynamical phase diagram (see the reviews 0,111 and references therein). The simplest model within this class 
of systems is the directed polymer in a two-dimensional random medium (see Q for a review), first introduced to 
model an interface in the low-temperature phase of two-dimensional disordered ferromagnets [10| . The statics of 
this model is rather well understood with exactly known critical exponents : at any temperature T, the directed 
polymer is in a disordered-dominated phase characterized by the roughness exponent £ eg — 2/3 and the droplet 
exponent 9 — 2( eq — 1 = 1/3 [9]. In comparison, the dynamics is much harder to study, and no exact result exists. 
Even numerically, the computational complexity changes completely between the statics that can be studied via exact 
transfer matrix methods that are polynomial in the size of the polymer 0] , and the dynamics where the determination 
of barriers is an NP-complete problem pd| . As a consequence, a complete characterization of the dynamical properties 
as a function of the temperature and external force (here we consider that the disorder strength is fixed) has remained 
a very challenging issue. An essential novelty with respect to the dynamics of the Sinai model described above, is the 
presence of the chain constraint of the polymer, so that even at zero temperature T = 0, the dynamics is already non 
trivial : the competition between interaction and disorder gives rise to a non-equilibrium depinning phase transition, 
between a pinned phase F < F^ ep where the interface remains blocked forever in some configuration, and a moving 
phase for F > F^ ep . This type of depinning transition has motivated a lot of experimental and theoretical studies 
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FIG. 1: Exactly known dynamical phase diagram as a function of the temperature T and the external force F of the one- 
dimensional Sinai model [2J : the important dimensionless parameter is a(F, T) = TF (the disorder strength has been fixed 
at some simple value). From the point of view of strong disorder renormalization [23 . I32I l33| . the logarithmic behavior 
x ~ ±(lni) 2 for F = corresponds to an "infinite disorder fixed point", whereas the anomalous diffusion phase x ~ t a for 
< a < 1 corresponds to a "strong disorder fixed point". A finite- velocity V > appears only for a > 1 , i.e. when the force 
is above some temperature-dependent value. 
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Two main questions: 

• Center-of-mass motion ? 

• Roughness at the largest scales ? 
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FIG. 2: The dynamical phase diagram for the directed polymer in a two-dimensional random medium as a function of the 
temperature T and the external force F is not exactly known. However one expects physically two very different regimes : 
(i) near the zero-temperature depinning critical point F^ ep , the transport at small temperature is dominated by depinning 
events, (ii) in the regime where T is finite and F is small, the roughness exponent ( eq = 2/3 of the equilibrium case F = is 
expected to describe the polymer roughness up to some large scale (that diverges as F — > 0). The transport is then dominated 
by thermal jumps of finite segments of the polymer between quasi-equilibrated metastable states. In this paper, we focus only 
on the regime (ii) and discuss the center-of-mass motion and the roughness at the largest scales. 
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(see the reviews [7|, |8[), and an exact mean-field solution has been obtained in [12j. The critical region with a small 
temperature around the depinning transition of an elastic line has been studied in [l3T ] . 

In the present paper, we focus on the opposite region of the dynamical phase diagram (see Fig. 2) where the 
temperature is finite and the external force is small. In this regime, one expects that the dynamics is dominated by 
thermal activations between locally quasi-equilibrated metastable states. In particular, the polymer is expected to 
keep during the motion its equilibrium roughness exponent ( eq = 2/3 up to some large scale (that diverges as F — > 0). 
This equilibrium roughness has been measured in dynamic simulations in 14], and in 15J when the temperature 
remains larger that some disorder-dependent threshold. However in other regimes, larger roughness exponents have 
been measured at larger scales [H, El- This issue of the roughness at the largest scales will be rediscussed later 
in the paper. But we first need to better understand the statistical properties of the barriers between the locally 
quasi-equilibrated metastable states. 



B. Usual creep scenario for finite temperature and small external force 

In the regime of finite temperature and small external force, the standard picture in the field (see the review [§| 
and references therein) seems to be that the dynamics then corresponds to a 'creep' motion, where the center-of-mass 
moves linearly in time as soon as the temperature is positive 

h G {t) ~ V(F,T)t (1) 
with a very non-linear force- velocity relation of the form 

K(T) 

V(F,T) ~ e-^ 1 (2) 

r — >U 

where the exponent fj, involves the dimension d of the interface (the directed polymer we are interested in corresponds 
to d = 1), the equilibrium roughness exponent £ eg , and the barrier exponent tp of the dynamics without external force 

(Ce, + d - 1p) 

Note that this relation is usually written in the form ^i usua i = (d — 2 + 2Q eq )/{2 — C, eq ) after using the additional 
assumption concerning the identity between the barrier exponent ip and the droplet exponent 9 = 2^ eq + d — 2. 
However, in our opinion, there is no convincing evidence of the equality ip = 9 neither theoretically (see [13] for a 
recent summary on the debate) nor numerically (the best numerical results presently available [l8l ] points towards an 
exponent of order ip ~ 0.49 rather different from the exactly known droplet exponent 9 — 1/3), we will consider in 
the following that the barrier ip is an independent exponent satisfying the bound 

tp>9 (4) 



C. Summary of the present paper : alternative scenario with a zero- velocity phase at small external force 

The usual qualitative argument given in favor of a finite velocity in Eq. [1] as soon as T > is that 'any barrier 
can be passed by thermal activation'. This is of course true, but this is not sufficient to conclude that the asymptotic 
velocity is finite, since the motion can also exhibit anomalous behavior as shown by the Sinai model at small external 
force (see the beginning of the introduction). To determine the asymptotic motion of the center-of mass at large 
times, one needs to study the time needed to travel over a given large distance. The more quantitative argument in 
favor of the finite velocity of Eq. [T] is a scaling argument where the barrier landscape existing in the absence of bias 
is 'tilted' to take into account the bias contribution : this scaling argument yields that the relevant barriers for the 
large-time dynamics corresponds to a barrier scale B*(F,T) and to a length scale l*(F,T) that are finite as soon as 
F > (and that diverge as F — * 0). We fully agree with this scaling argument that will be rediscussed below in section 
IIII CI However again, this is not sufficient to conclude that the velocity is finite, since exactly the same argument 
can be made for the Sinai model where an anomalous diffusion phase exists (see section IIII CI for more details) . The 
crucial property to determine whether the velocity is finite or not is the probability distribution of barriers B around 
this typical value B*(F,T). In the present paper, we explain that within the strong disorder renormalization in 
configuration space that we have introduced recently [l9| (see [2fj| for a more detailed presentation), one obtains an 
exponential tail for the distribution of renormalized barriers 



P(B) oc e B '^- T > 



(5) 



4 



As is well known in the field of disordered systems since Derrida's Random Energy Model I2ll . other disordered 
models sharing the same low-energy states statistics [13, and related Bouchaud's trap Models[2J] (see also [25[ for 
one of the first mention of exponentially distributed barriers in connection with extremal statistics), this 'innocent' 
exponential distribution for the barriers corresponds via the change of variable r = e B to a very broad power-law 
decay for the distribution of the trapping times r 

a(F,T) 

P ( T ) T « +OO T i+*(F,T) ( 6 ) 

The exponent depends continuously on the external force and on the temperature (we consider here that the disorder 
remains fixed) 

Since the characteristic barrier scale B* (F, T) grows with F and diverges as F — > as 

B*(F,T) cx — (8) 



the exponent a(F, T) vanishes as F — > as 



a(F,T) oc F» (9) 



As a consequence, the region of small external force where a(F, T) < 1 corresponds to a very broad distribution of 
trapping times with a diverging averaged value r = +oo. This invalidates the usual creep formula that is obtained by 
replacing all trapping times by the typical value T typ ~ e Btyp ~ e B ( F ' T ^ to obtain V ~ 1/ t *3/p (see Eq. [2]). We obtain 
instead that for a(F,T) < 1, the center-of-mass displacement then grows only sub-linearly in time 

h G (t) ~ t a(F ' T) (10) 

t — >oo 

i.e. the asymptotic velocity vanishes V = 0. 

Although the vast majority of papers on the subject never mentions the possibility of a zero-velocity phase, we are 
aware of three papers where the question of the probability distribution of barriers around the typical value has been 
raised : 

(i) twenty years ago, Ioffe and Vinokur [26] have proposed the anomalous sub-linear motion of Eq. [6] But for reasons 
that are very unclear to us, this possibility seems to have completely disappeared in the more recent literature. 

(ii) fifteen years ago, Bouchaud has pointed out in [27J that the usual creep argument might lead to incorrect results 
because of problem of strong fluctuations, and cites the Sinai model as an example where the transport is indeed 
harder than anticipated from the typical barrier alone. 

(ii) ten years ago, Vinokur, Marchetti and Chen [28| have proposed an exponential distribution of barriers, based 
on extremal statistics argument, and the corresponding power-law distribution for trapping times. However, these 
authors have supplemented this power-law distribution by a sharp cut-off to recover the usual finite velocity behavior 
of Eq. [T] In our opinion, there is no good reason to impose this sharp cut-off, since the same procedure for the exactly 
soluble Sinai model would give a wrong answer. 



D. Organization of the paper 



To summarize this long introduction, the aim of the present paper is to justify the existence of the anomalous 
diffusion phase of Eq. [TO] in the region of the phase diagram corresponding to finite temperature and small external 
force (see Fig. 2), via the use of some strong disorder renormalization procedure on the barriers. The paper is 
organized as follows. In section ITT] we briefly recall the main idea of strong disorder renormalization, and we describe 
the statistical properties of the "infinite disorder" fixed point that describes the thermal dynamics of the directed 
polymer without external force. In section 11111 we explain how the presence of a small external force F transforms this 
"infinite disorder fixed point" into a "strong disorder fixed point" with an exponential distribution of renormalized 
barriers, and we describe the consequences for distribution of trapping times and the large-time dynamics. In Section 
IIV1 we present detailed numerical simulations for a directed polymer driven in a traps landscape. In Section [Vj we 
compare our results concerning the roughness with previous works. Our conclusions are summarized in section IVII 
In Appendix A, we recall the subtleties associated to the type of the interactions (metric constraint versus elastic 
energy) to justify our choice to consider the metric constraint. 
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II. PROPERTIES OF THE 'INFINITE DISORDER FIXED POINT' FOR F = 

A. Strong disorder renormalization in configuration space 

Strong disorder renormalization (see for a review) is a very specific type of renormalization (RG) that has 
appeared in the field of quantum spin chains : this approach introduced by Ma and Dasgupta 30] has been developed 
by D.S. Fisher [3l[, who has introduced the crucial notion of "infinite disorder" fixed points where the method 
becomes asymptotically exact, and who has shown how to obtain explicit exact results for critical exponents and 
scaling functions for zero-temperature quantum critical points. This method has thus generated a lot of activity 
for various disordered quantum models [29$ . It has been then successfully applied to various classical disordered 
dynamical models, such as random walks in random media [321. [33|, reaction-diffusion in a random medium [3~jj . 
coarsening dynamics of classical spin chains ]35l , trap models [36j | , random vibrational networks [37| absorbing state 
phase transitions [H|, zero range processes [39j | and exclusion processes [4(| ■ In all these cases, the strong disorder 
RG rules have been formulated in real space, with specific rules depending on the problem. For more complex systems 
where the formulation of strong disorder RG rules has not been possible in real space, we have recently proposed in 
[T^ | a strong disorder RG procedure in configuration space that can be defined for any master equation as we now 
recall. 

The starting point of strong disorder renormalization in configuration space is the Master Equation describing the 
evolution of the probability Pt(C) to be in a configuration C at time t 

=J2 P t (O W ( C ' -> C ) - P t ( C ) W out (C) (11) 

at c , 

The notation W (C — > C) represents the transition rate from configuration C to C, and the notation 

W out (C) = J2 W (C - C) (12) 
c 

represents the total exit rate out of configuration C. 

For dynamical models, the aim of any renormalization procedure is to integrate over 'fast ' processes to obtain 
effective properties of 'slow' processes. The general idea of 'strong renormalization' for dynamical models consists in 
eliminating iteratively the 'fastest' process. The RG procedure introduced in [T^] consists in the iterative elimination of 
the state with the highest exit rate. We refer to (19LI20I] for a detailed description and derivation of the renormalization 
rules on the transition rates. Their most important property is their multiplicative structure that suggests that for a 
very broad class of disordered systems, the distribution of renormalized exit barriers 

B out = -\nW out (13) 

will become broader and broader upon iteration, so that the strong disorder renormalization procedure should become 
asymptotically exact at large time scales. Note that a very important advantage of this formulation in terms of the 
transition rates of the master equation is that the renormalized barriers take into account the true 'barriers' of the 
dynamics, whatever their origin which can be either energetic or entropic. 

After this general presentation (see [HI [2(| for more details), we now turn to the specific problem of a directed 
polymer in a two-dimensional random medium we are interested in. For the dynamics without external force, we have 
followed numerically the RG flow of the renormalized transition rates and we have found some "infinite disorder" 
fixed point (l9l . I20I ] as we now explain. 

B. Distribution of renormalized exit barriers at large scale 

The RG scale T is defined as the scale of the last eliminated exit barrier. At time t, the appropriate RG scale is 
thus 

r(i)=lnf (14) 

so that all metastable states of exit barrier B out > T, i.e. of exit time r = e Bout > t have been kept, whereas all 
metastable states of exit time r = e Bmlt < t have been eliminated. At large scale, one expects that the probability 
distribution of the remaining exit barriers B out > T will converge towards some scaling form 

Pr(B out ) ^ ' p(^i-L) (15) 

r^oo cr(r) v n r ) / 
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where P is the fixed point probability distribution, and where <r(r) is the appropriate scaling factor that represents 
the width of the renormalized distribution. 

For the directed polymer in a two-dimensional random medium, we have obtained numerically in (l9l . I20I ] that the 
width cr(r) grows asymptotically linearly with the RG scale V 

a(T) ~ T (16) 
r^oo 

and that the rescaled distribution of Eq. [T5] is extremely close to the exponential form 

P(x) ~ e~ x (17) 

Note that these two properties seem extremely robust within strong disorder RG since they hold for exactly in soluble 
models in d = 1 (29l | and have been also found numerically in quantum models in dimension d > 1 fill ]. 

C. Scaling between barriers and length scales 

To describe the thermal relaxation of disordered systems towards equilibrium when starting at time t = from a 
non-equilibrium initial state, it is useful to introduce the notion of some coherence length Wit) that grows slowly in 
time (see the reviews [42], |43| and references therein) . This coherence length separates the smaller lengths I < W (t) 
which are quasi-equilibrated at time t from the larger lengths I > Wit) which are completely out of equilibrium. 
Then full equilibrium is reached only when the coherence length reaches the macroscopic linear size Witeq) — L of 
the system. Within the droplet scaling theory proposed both for spin-glasses 0, SH and for directed polymers in 
random media [46] , the barriers grow as a power law of the length I 

B{l)~ft (18) 

with some barrier exponent ip > 0. The typical time t typ (L) associated to scale I grows as an exponential \nt ty p(l) ~ 
B(l) ~ ft. Equivalently, the coherence length-scale Wit) associated to time t grows only logarithmically in time 

l T {t) ~ (hit)* (19) 

In the numerical study of the relaxation towards equilibrium of an elastic chain in a two dimensional random medium, 
starting from a straight line at t = 0, this coherence length can be extracted from the behavior of the structure 
factor as a function of time [l|j], and the corresponding measure of the barrier exponent yields ip ~ 0.49 [l|| (see the 
comments on this exponent ip before Eq. [5]). 

Within the str ong disorder renormalization procedure, one may either extract a coherence length Zp as a function 
of the RG scale |2C| . or consider the statistics of the barriers associated to a given length [l9|, [2(| , and one obtains 
the scaling 

B(l)^A(T)ftu (20) 

where u is a random variable of order one. The barrier exponent tp obtained in [20j is of order ip ~ 0.47, i.e. close 
to the estimation tp ~ 0.49 measured in (l8j via Langevin dynamics. A systematic study of the dependence of the 
prefactor A(T) upon temperature is not yet available. This prefactor would be of the form C/T if the barriers were 
purely energetic, but since the strong disorder renormalization is defined on the transition rates, the prefactor A(T) 
may also contain entropic contributions and then have a more complicated temperature dependence. 

D. Final picture for the dynamics at F — 

The final picture for the dynamics of the directed polymer of length L at finite temperature T and no external force 
F = is thus the following : 

(a) during the regime ?t(X) < L where the coherence length lr(t) is smaller than the length L of the polymer, the 
polymer is quasi-equilibrated with its equilibrium roughness exponent £ eq = 2/3 only on smaller length scale than the 
coherence length I < lr{t), whereas larger length scale I > Wit) are still completely out-of-equilibrium. The "infinite 
disorder fixed point" describes the hierarchical structure of growing barriers that have to be passed to equilibrate on 
larger and larger length scales, and this is why the coherence length grows only logarithmically in time 

™~{mf (21) 
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(b) when the coherence length lr(t) reaches the length L of the polymer, the full polymer is characterized by the 
equilibrium roughness exponent C, eq — 2/3. The coherence length lx{t) cannot grow anymore, so that the renormal- 
ization procedure has to be stopped at the RG scale of order 

T L ~ AiT)^ (22) 

One then expects that the polymer of roughness exponent Q eq = 2/3 will become able to move slowly in the transversal 
direction between metastable states that are separated by distances of order i> e « and by barriers distributed with the 
fixed point distribution ( Eq. [T51 [TH] and [T7|) for the scale V l , that presents the exponential decay 

Pv L (Bout) K - =U"^f (23) 

On the renormalized scale Tl, one expects that the transverse motion of the center-of-mass corresponds to an effective 
one-dimensional Sinai model, where the unit distance scale is of order L^" q and the unit barrier scale is Yl ~ A(T)L^. 
As a consequence for lni 3> A(T)L^, one expects the logarithmically slow behavior 

M*)=.±i«-( 3 ^)' m 



III. PROPERTIES OF THE 'STRONG DISORDER FIXED POINT' FOR SMALL FORCE F 

A. Analysis of the RG flow for small F 

The "infinite disorder fixed point" described above for F = characterizes some 'criticality' in the time direction in 
the following sense : there is no characteristic scale for the barriers except the RG scale T itself. This scale invariance 
will be broken by the introduction of some external force F that will introduce some characteristic length-scale and 
thus some corresponding scale B*(F, T) for the barriers. However, if the external force is very small, the scale B*(F, T) 
(that diverges as F — > 0) will be very large. We thus expect that the RG flow can be analyzed in terms of two regimes 

(i) during the first regime 1 <C T < B* (F, T) , the distribution of the renormalized exit barriers will follow the same 
RG flow as in the absence of force, i.e. it will converge towards the scaling of Eq. [15l with the same exponential 
rescaled distribution of Eq. [17] 

Pr,F, T (Bout) * .(r^) 6 "™ ( 25 ) 

with a width a(T, F, T) which grows linearly in T as in Eq. |T5] 

aCT.F.T) ~ T (26) 
i«r«B«(f,T) 

(ii) when the large scale T ~ B*(F,T) is reached, the width saturates at the finite large value B*(F,T) 

a(T,F,T) ~ B*(F,T) (27) 

r>B*(_F,T) 

instead of the flow towards infinity that characterizes the critical case F = 0. Since this width remains finite asymp- 
totically as r — > oo, one speaks of a 'finite-disorder fixed point'. However, since this width B*(F,T) diverges at small 
force , the region of small F is a 'strong disorder fixed point', where the asymptotic accuracy of the renormalization 
approach is of order 1/ B*(F, T). This notion is thus very useful to study the vicinity of 'infinite disorder fixed point' 
in the space of parameters and we refer to the review [29] for more detailed discussions. In one-dimension, where the 
strong disorder RG procedure can be followed exactly, the crossover of the width between the regimes (i) and (ii) 

described above, is of the form (translated in our present notation) : a(T,F,T) = B*{F 1 T) 1 - e~ B *(^. T > [H,[3l|. 

For our present analysis, the important point is that the saturation scale B*(F,T) is large enough, so that the RG 
flow during the first regime 1 C T « B*(F,T), that behaves as the critical flow, contains sufficiently RG steps to 
have converged towards the scaling form of Eq. [25l So, when saturation occurs at scale B*(F,T), the probability 
distribution of renormalized exit barrier follows the exponential form 

Pr{B) ^ B*(F,T) (28) 
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B. Physical meaning of the saturation 



The physical meaning of the two regimes described above is as follows. During the hrst regime (i), the external 
force F is so small that the flow is very similar to the flow for F = 0. In particular, the barriers grows upon iteration 
on scales 1 <C T <C B*(F,T) and the motion is not yet directed along the bias. When the saturation occurs at the 
large scale B*(F, T), this means on the contrary that the motion becomes effectively directed in the direction of the 
external force F for scales T > B*(F,T). The barriers against the bias are not renormalized anymore, and one can 
stop the renormalization procedure. The appropriate model on this scale is then a directed model along the bias, 
with barriers distributed as in Eq. [221 As recalled in the introduction (see the discussion between Eqs. [5]and[()]), the 
corresponding distribution of the trapping time t = e B is then a broad power law 

P (r) - ^W7F) ( 29 ) 

with exponent 

a(F,T) = (30) 



C. Determination of the saturation scale B*(F,T) 

Within our present RG framework, the saturation scale B*(F, T) should be determined as the limiting value of the 
width (Eq. I27p of the renormalized distribution of barriers. To determine the dependence on the external force at 
small F, we may rephrase the scaling argument which is usually used in the field |8j as follows. On the length scale 
the barriers in the absence of external force follow the scaling of Eq. 1201 with the barrier exponent "0. If the external 
force F is small, one may take into account its effects by a 'tilt' of the landscape that lowers the barriers against the 
force in the following way 

Bit) ~ A(T)l^u- (31) 

where £ er? represents the equilibrium roughness exponent, so that the correction in Eq. [31] corresponds to a transversal 
move of order for the segment of length I of the polymer. Since (1 + Q eq ) > ip, the force term always dominates 
at sufficiently large length scale. The length scale l*(F, T) that will give rise to the biggest barriers can be obtained 
by differentiating Eq. [31] with respect to I. Dropping constants of order 0(1) one obtains the length scale 

/TA(T)\ (e<s,+i-« 



and the corresponding barrier scale 



B*(F,T) F ^ Q A(T)(^iy (33) 



where (i is the exponent that usually appear in the creep formula (see Eq. [3] and associated comments) 

" ~ (C eq + l - V) 

So the power-law exponent of Eq. [30] vanishes for F — » as 



(34) 



i i / f y 

In the case of the Sinai model, one may actually use exactly the same arguments with the following changes : 
the length / is now along the direction of the motion, the barrier exponent -0 = 1/2 simply describes the Brownian 
fluctuation of the random energy landscape, the prefactor A(T) has the simple T-dependence A(T) = C/T because 
the barriers are purely energetic, and there is no roughness exponent Ceq = since the model concerns a single particle. 
For this special case, one obtains \i — 1 and the following expression for the power-law exponent 

(F, T) oc TF (36) 



in agreement with the exact results mentioned at the beginning of the introduction. 
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D. Final picture for the dynamics at small external force F 



In the case of the Sinai model, the saturation scales l*(F,T) for the length and B*(F,T) for the barrier are the 
only finite scales present in the problem. The strong disorder renormalization procedure should then be stopped 
at the scale B* (F, T) , and one ends up with an effective directed trap model. We refer to [33| for a more detailed 
presentation of the quantitative relations that can be derived at large scales between the Sinai model in an external 
force and the one-dimensional directed trap model. 

For the present model concerning the directed polymer, the polymer length L introduces another length scale in 
the problem. Of course, one is interested into large polymer length L»l, but since the length l*(F,T) diverges as 
F — ► (Eq. [32]) . one needs to distinguish various regimes in terms of the polymer length L 



1. Regime 1< L < 1*{F,T) 



In the regime lci« l*(F,T), the relevant barriers for the motion of the polymer as a whole corresponds to the 
RG scale of Eq. [22] which is well below the saturation scale B*(F,T) introduced by the force F. This means 
that the force F is then too small to impose a directed motion between two neighboring quasi-equilibrated metastable 
configurations. The motion will become effectively directed only on larger transversal lengthscales separating many 
metastable configurations. 



2. Regime L ~ l*(F,T) 



For a polymer of length of order L ~ l*(F, T), the motion of the polymer corresponds to an effectively directed trap 
model. Each trap of trapping time Tj corresponding to a transverse displacements of order Xi — (l*(F, T)pXi (where 
Xi is a random variable of order 0(1) with a finite averaged value Xi < oo). After n traps, the total displacement of 
the center of mass will be of order 

h G (n) = Xl +x 2 + ... + x n ~ n[l*(F,T)] C (37) 

n — >oo 

whereas the total time needed to escape from these n traps will be of order 

t(n) = ti + t 2 + ... + t„ (38) 

As is well known in the field of Levy statistics [2| , the sum of n variables distributed with the broad power-law of Eq. 
[29] grows linearly in n only when the average value f remains finite, i.e. for a > 1 

tin) ~ rif if a(F, T) > 1 (39) 

n — >oo 

However for a < 1, the average value of the trapping time diverges r = +oo and the sum grows more rapidly as 

t{n) ~ n« if a(F, T) < 1 (40) 

n — >oo 

It is then clear that the asymptotic velocity is finite only for a(F,T) > 1. For a(F,T) < I, the velocity vanishes 
V = and the center-of-mass displacement grows only sub-linearly in time 

h G (t) ~ t a{F -^ if a(F,T) < 1 (41) 

t— »OG 



3. Regime L> V (F,T) 



Finally in the regime L 3> l*(F, T), the strong disorder RG procedure has to be stopped at the RG scale B*(F, T) 
of Eq. l33l corresponding to the length scale l*(F, T) of Eq. [32] The picture that emerges is then some kind of directed 
'parallel' trap model, where a number of order L/l*(F, T) segments of typical length l*(F, T) have to move in parallel 
in a trap landscape where the statistics of trapping times is given by the broad power-law of Eq. [29] with exponent 
given in Eq. 1351 The main question is then to determine the form of the effective interaction between these large 
segments from the knowledge of the microscopic interactions between monomers. This issue is discussed in detail 
in Appendix fA] where we first recall why the metric constraint is a safe choice for the microscopic interactions (in 



10 



contrast with a quadratic elastic energy that is known to lead to an unphysical roughness exponent in some cases), and 
where we then explain why the metric constraint is also an appropriate choice for the effective interactions between 
consecutive segments of size l*(F, T). 
From the two properties : 

(i) each segment of size I* (F, T) would follow the sub- linear motion of Eq. [4J] if it were alone 

(ii) the metric constraint between these segments can only delay (but not accelerate) each segment with respect to 
the motion it would have followed if it had been alone 

we conclude that the center-of-mass motion will also follow the asymptotic sublinear behavior of Eq. [?T] 

h G (t) ~ t a ^ if a{F,T) < 1 (42) 

t — >oo 

This is confirmed by the numerical simulations presented in the next section, where we moreover discuss the behavior 
of the interface width at the largest scales (i.e. at scales much larger than l*(F,T) for the initial microscopic model). 



IV. NUMERICAL STUDY OF A DIRECTED POLYMER DRIVEN IN A TRAPS LANDSCAPE 

In this section, we define and study an effective directed model within a landscape made of traps, which is expected 
to be the appropriate coarse-grained model of the true microscopic non-directed dynamics in the regime L » l*(F,T) 
discussed above in section IIIID 31 We stress that a single monomer of this effective model represents a quasi- 
equilibrated segment of length I* (F, T) of internal roughness Q eq — 2/3 of the true microscopic model (see the previous 
sections for more explanations). 



A. Definition of the effective directed model within a traps landscape 

We consider a directed polymer of length L defined by the heights (hi, h 2 , ■ hi) with cyclic boundary conditions 
hL+i = h\ and with the metric constraint 

\h i+ i-hi\ = ±l (43) 
The initial configuration is the 'flat' zig-zag configuration 

hn = 1 

h 2i+ i = (44) 

The traps landscape is defined as follows : the random trapping times r(i,hi) are independent and drawn from the 
power-law distribution 



q (r)=9(r>l)-^ (45) 



The dynamics is defined by a Master Equation of the form of Eq. [TT] To respect the chain constraints of Eq. l43l 
the monomer (i) of the configuration C = (hi, hi,..., hi)) is 'movable' only if hi—\ = hi + 1 = /ij+i, and the possible 
movement is then hi — ► hi + 2. The total rate out of the configuration C is thus given by 

L 

W out (C = {hi,h 2 ,...,h L }) =J2Wt(C) (46) 

i=l 

W l (C = {hi,h 2 ,...,h L }) = 6h *- h * +1 (47) 

T[i,hi) 

The two main observables in the field of interface dynamics arc 49] 

(i) the height hc(t; L) of the center-of-mass of an interface of length L as a function of time t 



1 L 

h G (t;L) = I J2h i (t) (48) 



i=l 



(ii) the width w(t) of the interface of length L as a function of time t 



■• 2 (t;L) = j^2[h i (t)-h G (t)} 2 (49) 
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B. Numerical details 



In our numerical study, we have used the 'Bortz-Kalos-Lebowitz algorithm' [47| which is a faster-than-the-clock 
algorithm where each program iteration computes the time and the site where the next movement occurs [48j . The 
idea is that from the knowledge of W out , the escape time t esc from configuration C is a random variable drawn from 
the law 



P£ xlt (t e 







Then the monomer (i) which is effectively moved hi 

?r(i) 



W oat (C)e- w ^ c ^ 

-» hi + 2 is drawn with the probability 

W t (C) 
W out {C) 



(50) 



(51) 



™ ax wished at the end of the simulation fixes the 



In this framework, the total displacement he of the center-of-mass is the natural variable, and one computes for a 
given dynamics in a given disordered sample the time t(ha) needed to reach a given displacement ha of the center- 
of-mass. Similarly, one computes the width w(hG) of the interface for a given displacement he of the center-of-mass 
starting from the flat initial condition of Eq. 1441 The numerical results presented below correspond to measures at 
the values ho — 2, 4, ....K^ ax : where the maximal displacement h 
CPU time. 

For the case a — 0.5 where the averaged trapping time diverges r = +oo ( Eq l4"ij)). we have made studies with 
fomax ^ anc j w ith hQ ax ~ L 15 for the following sizes and the corresponding number n s (L) of disordered samples : 

(i) for h~Q ax ~ L, we have studied sizes in the range 10 3 < L < 24. 10 3 with a statistics between n s (L = 10 3 ) = 4.10 5 
and n s (L = 24.10 3 ) = 1600. 

(ii) for Kq clx ~ L 15 , we have studied sizes in the range 10 2 < L < 4.10 3 with a statistics between n s (L — 10 2 ) = 
3.10 6 and n s (L = 4.10 3 ) = 500. 

For the case a = 2 where the averaged trapping time is finite r < +oo ( Eq I45[) . we have made studies with 
h™ ax ~ L for sizes between 10 2 < L < 10 4 with a statistics between n s (L = 10 2 ) = 4.10 5 and n s {L = 10 4 ) = 1600. 
(Note that all other parameters being the same, the CPU time turns out to be larger for the case /x = 2 than for the 
case fi = 0.5, because the number of 'movable monomers' in the stationary regime is larger for the case fi = 2) 

We give below results for single histories in a given disordere d sample, as well as disorder-averaged results denoted 
by an overbar. In particular, we will be interested in Int(hc) that represents the disorder-average of the logarithm 
of the time needed to reach a displacement he of the center-of-mass. For the width defined in Eq. 221 we moreover 
introduce the following notation 



w av (t,L) = (w 2 (t;L) 



1/2 



(52) 
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FIG. 3: (Color online) Typical dynamics of a single interface : log- log plot of the center-of-mass height ha(t) and of the 
interface width w(t) as a function of the time t time for a polymer of length L = 2000 (a) for the case a = 0.5 where the 
averaged trapping time diverges f = +00 ( Eq |45fl . the dynamics is very intermittent, (b) for the case a — 2 where the averaged 
trapping time is finite f < +00 ( Eq I45[) . the center-of-mass motion is smooth. 
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C. Typical dynamics of a single interface 



We first consider a single history in a given disordered sample for a polymer of length L = 2000. The dynamics of 
the center-of-mass position ha{t) (Eq. and of the width w(t) (Eq. |4T))) are shown on Fig. [31(a) for the case a = 0.5 
where the averaged trapping time diverges r = +oo ( Eq [45|) . and on Fig. [3](b) for the case a — 2 where the averaged 
trapping time is finite r < +oo ( Eq [45)l . This comparison shows that the dynamics is very different in the two cases. 
For a = 0.5, the stationary regime is characterized by a very intermittent dynamics of the center-of-mass hc(t) and 
of the width w(i), which remain pinned for long time intervals separated by rapid avalanches. On the contrary for 
a = 2, the motion of the center-of-mass is smooth, and in the stationary regime, the width w(t) fluctuates rapidly 
around its time-averaged value. 



h(i) 









(a) 



Mi) 



500 



1000 



1500 



2000 




2000 



FIG. 4: (Color online) Typical configurations in the stationary regime for a polymer of length L — 2000 (a) for the case a — 0.5 
where the averaged trapping time diverges r = +oo ( Eq I45[) . very large regions of slope 1 are present (b) for the case a = 2 
where the averaged trapping time is finite r < +oo ( Eq I45[) , there are more structures on shorter scales. 



We show on Fig. 0] typical corresponding configurations in the stationary regime. For a — 0.5, one clearly see on 
Fig. |4] (a) very large regions of slope 1, whereas for a — 2, the configuration of Fig. [4] (b) presents more structures on 
shorter scales. 



D. Disorder-averaged behaviors for a given size L 

We have shown above on Fig[3]the dynamics of hc(t) and w(t) for a single history of a polymer of size L. We show 
on Fig. [5] the same observables after averaging over n s = 2000 disordered samples. 

We note t*(L) the crossover time between the transient regime where the disorder-averaged width u> al ,(t, L) ( Eq. 
I52[) grows and the stationary regime where the disorder-averaged width saturates towards a time-independent value 

w av (t,L) ~ w sat {L) (53) 

t»t*(L) 

For a — 0.5, we find that in the stationary regime, the plot of In ha as a function of \nt(ho) (see Fig. [5] a) 
corresponds to a slope 0.5 

lnh G ~ 0.51ni+ ... = aln< + ... (54) 
t>t*(L) 

in agreement with the sub-linear motion predicted in Eq. 1421 Of course for a = 2 we find that in the stationary 
regime, the plot of \nhc as a function of hit (see Fig. [5]b) corresponds to a slope 1 corresponding to the usual case 
of finite velocity hQ(t) ~ Vt. 

For a — 0.5, we moreover measure that, in the initial transient regime, the center-of-mass and the width both grows 
with the same exponent [3 in time 

lnh G (t,L) ~ /31nt+... 

t«t*(i) 

]nw av (t,L) ~ 0lnt+... (55) 
«t»(i) 
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FIG. 5: (Color online) Disorder-averaged properties of the center-of-mass height ha and of the interface width w as a function 
of time (log-log plot) : data averaged over n s = 2000 disordered samples for a polymer of length L = 2000 (a) for the case 
a = 0.5 where the averaged trapping time diverges f — +oo ( Eq I45p . the center-of-mass motion is sublinear with exponent 
(3 ~ 0.25 in the transient regime (Eq. I55[) and with exponent a = 0.5 in the stationary regime (Eq. I54|l (b) comparison with 
the case a = 2 where the averaged trapping time is finite r < +oo ( Eq I45|l . 



with a value of order 



that will be interpreted in section HV Gl below. 



(3 ~ 0.24 



(56) 



E. Roughness exponent in the stationary regime 




FIG. 6: (Color online) Dynamics of the width of the interface for a — 0.5 : (a) lnui a „ as a function of ln/iG for various sizes 
L (b) In Wsat (L) as a function of InL : the slope corresponds to the roughness exponent f = 1. 



We now consider the dynamics of the disorder-averaged width defined in Eq. [52] as a function of the center-of-mass 
displacement \iq. 

For a = 0.5, we show on Fig. [6] (a) the dynamics of the width w av for various polymer sizes L. We show on Fig. [6] 
(b) the log-log plot of the saturation value w sat {L) defined in Eq. as a function of L : we find a slope of order 1 

\nw sat {L) ~lnL + ... (57) 

Our conclusion is thus that for a < 1 , the roughness exponent is £ = 1 , as already suggested by the shape of typical 
configurations in the stationary regime (see Fig. S]a ) 
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For comparison, we show on Fig. [7]the same observables for the case a — 2. In particular, the log-log plot of Fig. [7] 
(b) corresponds to a slope of order In W sa t(L) ~ 0.95 In L for the width in the stationary regime. This high value can 
be understood from the typical configuration in the stationary regime shown on Fig. |4]b. This suggests that even in 
the finite velocity phase a > 1 where the averaged trapping time is finite (Eq I45[) . it is very difficult for the different 
regions of the polymer to remain synchronized during the dynamics, so that the the polymer tends to be stretched on 
large distances with a slope close to the maximal slope 1 allowed in the presence of the metric constraint. 



F. Finite-size scaling forms in t and L in the case a < 1 




FIG. 8: (Color online) Finite-size scaling analysis in t and L the width of the interface in the case a = 0.5 : (a) lnw av as 
a function of Int for various sizes L (b) Data collapse of the same data using the variables In (w sa t{L)/L) as a function of 
In (t/L z ) (see Eq. EHJ with 2 = 4. 



In the field of interface dynamics, it is usual to summarize the crossover between the transient dynamics and the 
stationary dynamics by the following scaling form for the width [49] 

w av (t-,L)~L<<p(jA (58) 

where £ represents the roughness in the stationary regime, where z is the dynamical exponent of the crossover time 
t*(L) ~ L z and where the scaling function <f>(x) has the following behaviors : 



15 



(i) it converges towards a constant at infinity (j){x — > oo) ~ const to recover the saturation value (Eq. I53() at large 
time ; 

(ii) it behaves as a power-law at small argument : <p(x — > 0) ~ x 13 where the exponent j3 governs the growth in time 
during the initial transient behavior 

w av (t;L) ~ tf-P'tPrvt 13 (59) 

where the scaling relation Q — [3z between exponents is expected to obtain a L-independent behavior in the transcient 
regime (49| . 

We show on Fig. [5] that a good data collapse can be achieved with the values £ = 1 and 

z(a = 0.5) ~ 4 (60) 

The corresponding exponent of the scaling function <fi(x) at small argument is measured to be of order 

0{a = 0.5) ~ 0.24 (61) 

in agreement with the previous estimate of Eq. 1561 

Equivalently, as a function of the center-of-mass position h G , the width satisfies the finite-size scaling form 

w av {h G \L)~Li>(^\ (62) 

where the scaling function ip(x) has for asymptotic behaviors ip(x — ^ oo) = const and tp(x — > 0) oc x, so that in the 
transient regime, one has the simple proportionality 

w av (h G ;L) ~ (const) h G (63) 

The compatibility between the finite-size scaling forms of Eq. [58] and of Eq. [62] yields the following finite-size 
scaling form for the center-of-mass position 

h G {t;L)~LH(±-\ (64) 



L z 



where H(x) is a scaling function (data not shown). 



G. Simple arguments to determine the exponents for < a < 1 

We expect that the roughness C = 1 found above numerically for the special case a = 1/2 remains the same in the 
whole phase < a < 1 where strong non-self-averaging effects of the trapping times (Eq |4"5|) occur 

C(0 < a < 1) = 1 (65) 

We may now use a simple argument to determine the other exponents as a function of a in the interval < a < 1. 
For h G large, the time t(ha,L) needed to make a center-of-mass displacement of order h G for a polymer of length 
L may be estimated from the maximal T max (N) trapping time among N — (Lh G ) independent variables drawn with 
the probability distribution of Eq [431 

Tmax(N) -N 1 ^ ~ (LM 17 " (66) 

From this maximal trapping time, the minimal rate transition rate encountered is from Eq. 1471 

W m T (671 

Assuming that the time t scales as the inverse of this minimal rate t ~ 1/W™ l t n ~ x (Lha) 1 / 01 , we obtain by inversion 
the following scaling in the stationary regime 

h G (t;L) ~ (68) 
t»t*(L) L 
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The scaling function H(x) of Eq. [64]thus grows asymptotically as the power-law H(x — > oo) ~ x a and the dynamical 
exponent z describing the scaling of the crossover time t*(L) ~ L z reads 

z(a) = - (69) 
a 

For a = 1/2 used in our simulations, this corresponds to z — 4 in agreement with the value measured above (EaslSOp. 

To determine the exponent f3 of the initial transient regime, the above argument has to be changed as follows. We 
expect that in the transient regime, an extensive number (pL) of monomers (with a finite density p) are still pinned 
in their initial positions. Equivalently, there is a number pL of segments of typical length l p ~ L/(pL) ~ 1/p, that 
have made a move of typical amplitude hi ~ I, because the dynamics is building the roughness £ = 1 at these small 
scales. As a consequence, the center-of-mass displacement is of order ho ~ pip ~ l p - Since all these segments are 
independent, we now have to evaluate the time ti p for a single segment of size l p to move over a distance of order l p , 
i.e. the time ti p to overcome l p trapping times. The maximal trapping time among these scales as (see Eq. I66[) 

W(tf p = £)~^ 1/0, ~ip /a (70) 

The time can be thus estimated as in Eq. [571 : i; p ~ l/W™J n ~ xl^ a . Using ftc ~ Z p derived above, we obtain by 
inversion the transient behavior 

h G (t;L) ~ (71) 
*<**(£) 

The scaling function -ff(x) of Eq. [64] thus grows at the origin as the power-law H(x — > 0) ~ x a l 2 i.e. the transient 
exponent f3 reads 

(3(a) = | (72) 

For a — 0.5, this is in agreement with the value (3 ~ 0.24 measured above (see Eqs IS"6l and l6"Tj) . Finally, with the 
values of Eqs l6"5l IBTH and [721 the transient behavior of the width of Eq. [591 simplifies into the same scaling behavior as 
the center-of-mass displacement ho (Eq. [7T|) 

w av (t;L) ~ t* (73) 

t<t*(L) 

as expected from Eq. [63j 

In conclusion, the dynamics in the zero-velocity phase < a < 1 can be understood in terms of simple arguments 
based on the statistics of trapping times. The critical exponents that are computed from these arguments are in 
agreement with our numerical simulations both in the stationary regime and the transient regime. 

V. COMPARISON WITH ROUGHNESS EXPONENTS MEASURED IN PREVIOUS WORKS 

For a polymer with the metric constraint at finite temperature driven by a small force, we may summarize the 
behavior of the roughness discussed in the above sections as follows : 

(i) at scales smaller than l*(T, F), the force is too small to really change the energy landscape seen by the polymer 
at F = 0, and thus the roughness is expected to be governed by the equilibrium value ( eq = 2/3. 

(ii) at scales larger than l*(T,F), we have found in our numerical study of the previous section that, at small force 
where the exponent of the trapping time power law (Eqs |2"91 and [3U)) is in the interval < a(F,T) < 1, the roughness 
becomes £ = 1. 

In this section, we compare this scenario with previous works. 

A. Relation with other works finding a roughness exponent £ = 1 in the presence of a metric constraint 

We are aware of two works where a roughness exponent C = 1 nas been measured for the dynamics in the presence 
of a metric constraint: 

(i) Sneppen has found a roughness exponent £ = 1 for the following self-organized interface dynamics called 'model 
A' in [50j : at each iteration, the site with the smallest pinning force among the sites that are allowed to move without 
breaking the metric constraint, moves forward by one unit and a new pinning force is drawn for this site. 
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Self-organized-criticality models are usually not directly related to finite-temperature dynamics. Here for instance, 
the important differences with the effective directed traps model defined by the master equation of Eq. [TT] with Eq. 
|47] are the following : in Sneppen's model, there is no real time, but only a number of iterations that corresponds to 
the total displacement Hq of the center-of-mass; and at each iteration, it is the always the least pinned movable site 
that moves, whereas in the master equation of Eq 111! the time and the site of the next moves are drawn with Eqs l50l 
and I5ip . Nevertheless, there is clearly a close relationship between the two models, in particular in the mechanism 
that generates a roughness exponent £ = 1 in the stationary state (see the typical configurations on Fig. 2] (a) of the 
present paper and Fig. la of Ref. |50j). 

(ii) Tang and Leschhorn have found in their analysis of the dynamics slightly above the directed percolation 
depinning threshold [5ll . |52| that the moving interface is not sclf-affinc, but is a mixture of two kinds of behaviors : 
there exists finite pinned segments that have for roughness the exponent of the depinning transition Q^ ep — Cdp ^ 0.63 
whereas segments with a slope of order 1 are found between them. 

Since the region discussed by Tang and Leschhorn is on the horizontal axis of Fig. 2 slightly above the depinning 
critical point, it is very far from the region near the vertical axis of Fig. 2 discussed in the present paper. However, 
the two pictures that emerge have nevertheless some similarities : in our case, the local roughness is determined by 
the equilibrium exponent Q eq = 2/3 up to the length /*, in their case the local roughness is the depinning roughness 
Cdep = (dp ~ 0.63, but in both cases the roughness becomes equal to 1 at largest scales, i.e. the maximal roughness 
which is physically acceptable. This seems to indicate that in both cases, it is very difficult for the different regions 
of the polymer to remain synchronized during the dynamics, so that the maximal slope 1 allowed by the metric 
constraint is actually reached to maintain the full polymer together. 

B. Roughness exponents measured in other works in the presence of an elastic energy 

Most of the numerical studies on the dynamics at finite temperature have been done for a Langevin dynamics in 
the presence of an elastic quadratic energy between monomers : 

(i) the equilibrium roughness exponent ( eq — 2/3 has been measured during the driven dynamics in [Til ], and in 

[la ] when the temperature remains larger that some disorder-dependent threshold. In our opinion, this means that 
in these two simulations, the l eng th L of the polymer was smaller or of the order of the length l*(F,T). 

(ii) at lower temperature in [15| , the authors have measured an effective roughness exponent of order £ ~ 0.9. In the 
study 16] concerning the limit T — > 0, the authors have found a crossover between the equilibrium value ( eq = 2/3 at 
small scales and the critical depinning value (d e p ~ 1-26 at larger scales, and have related this crossover to Functional 
RG calculations [53j . Since any roughness exponent £ > 1 is unphysical (see [54l |55| and the discussion in Appendix 
A) , our opinion is that it would be very interesting if the simulations of [161 ] based on a quadratic elastic energy were 
done with the metric constraint instead, to see what physically acceptable roughness exponent Q < 1 would actually 
emerge at large scales. 

VI. CONCLUSIONS AND PERSPECTIVES 

In this paper, we have considered the dynamics of the directed polymer in a two-dimensional random potential in 
the regime where the temperature T is finite and the external force F is small. We have explained how the "infinite 
disorder fixed point" that describes the dynamics for F — becomes a "strong disorder fixed point" for small F 
with an exponential distribution of renormalized barriers. Since the corresponding distribution of trapping times then 
only decays as a power-law P(t) ~ l/r 1+Q , where the exponent a(F,T) vanishes as a(F,T) oc F^ as F — » 0, we 
have concluded that the motion is only sub-linearly in time ha(t) ~ t a ( F,T > in the region a(F,T) < 1, i.e. that the 
asymptotic velocity vanishes V = 0, in contrast with the usual creep scenario where the velocity is finite as soon as 
(T > 0,F > 0). All along the paper, we have discussed the similarities with the Sinai model with bias, where an 
analogous zero- velocity phase has been established long ago by rigorous methods 0, 0, HI, @] and where the asymptotic 
exactness of the strong disorder renormalization has been demonstrated explicitly by a direct comparison with the 
available rigorous results [2{| H2, [Hj]. We have then checked the presence of the predicted zero- velocity phase by 
numerical simulations of a directed polymer with a metric constraint driven in a traps landscape. We have moreover 
obtained that the roughness, which is governed by the equilibrium exponent £ eg = 2/3 up to the large scale I*, is 
equal to £ = 1 at the largest scales. 

An important issue is of course whether such a zero- velocity phase also exists for interfaces of higher dimensionalities 
d > 1 driven in a random medium of dimension (d+1), and more generally for other classes of driven extended systems. 
Since "collective transport in random media is an impossibly broad subject" 0, a general answer clearly goes beyond 
the present work. However we think that the essential property needed to have a zero-velocity phase at small force 
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is the presence of a positive barrier exponent ip > for the dynamics at F = 0, and this should be the case for a 
broad class of disordered systems in finite dimensions. Then the dynamics for F — will be logarithmically slow and 
should correspond to some "infinite disorder fixed point" for the renormalized barriers, that transforms into a "strong 
disorder fixed point" at small F, with an exponentially distribution of renormalized barriers. Another argument in 
favor of this general scenario is the extremal statistics argument on the barriers [H|,[28| that also lead to an exponential 
tail for the probability distribution of large barriers, and thus to a power-law decay for the trapping time distribution. 
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APPENDIX A: TYPE OF INTERACTIONS : METRIC CONSTRAINT OR ELASTIC ENERGY 

In the text, the interaction used to maintain the continuity of the interface is a metric constraint, both for the 
microscopic model and for the effective directed trap model described in section UlID 31 In the present appendix, we 
explain this choice. 



1. Subtleties in the choice of microscopic interactions 



Let us first recall the two types of microscopic interactions that are usually considered in the literature for a polymer 
described by the heights {hi} of monomers : 



a. Metric constraint between the heights of two neighboring monomers : \h i+ i — h t \ < 1. 

This metric constraint is for instance very much used in numerical studies using transfer matrix methods (see the 
review Q). It can can be moreover justified in various microscopic models, in particular when the directed polymer 
represents an interface in a two-dimensional disordered ferromagnet well below the critical temperature T c , which was 
the original motivation to introduce the directed polymer model [Iol |. 



b. Elastic energy of the form : E e i(hi+i — hi) ~ — hi) n , with n — 2 usually. 



It is of course very common in physics to replace a hard constraint by a soft constraint. In the present case, the 
use of an elastic quadratic energy usually comes from a small-gradient expansion, an hypothesis which has to be 
consistent with the results obtained for the roughness. Indeed, the roughness exponent £ has to satisfy the bound 
( < 1 for at least two reasons (EH, : first, if one obtains C > 1, the gradient is not small, and thus there is an 
inconsistency with the small gradient expansion used to obtain the elastic quadratic energy ; second, if one obtains 
C > 1, the elastic energy diverges in the thermodynamic limit L — > oo. 



c. Discussion 



The metric constraint can be seen as the limit n — > oo of the power n of the elastic energy. As a consequence, there 
is room for various universality classes between n = 2 and n — oo. So one should not assume a priori that all these 
interactions are equivalent, but test for each case of interest if they lead to the same results or not. 

For the directed polymer at equilibrium at any temperature T and no external force in dimension 1 + 1, the same 
equilibrium roughness exponent C, eq = 2/3 arises if one considers the metric constraint or the elastic energy with n = 2 
(see the review[9(). However, at the zero-temperature depinning transition (see Fig. 2), some subtleties arise (see the 
detailed discussion in (5o| ) : it turns out that the model defined in terms of an elastic quadratic energy (case n = 2) 
gives a roughness exponent ( n =2 ~ 1-25 which is unphysical because it is strictly greater than one ( n=2 > 1 [EH. |55|. 
This is in contrast with the model defined either with the metric constraint or with an elastic energy of at least quartic 
order (n > 4) where the roughness exponent is physical ^metric _ 1 and actually takes a value ^metric ~ 0.63 ~ Cn=4 
[5o| that can be understood in terms of pinning by Directed Percolation clusters fo"ll |56| . 
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In conclusion, since the microscopic model with an elastic quadratic energy (case n = 2) is known to lead to 
unphysical results at the depinning transition (Cn=2 > 1), it is clearly not a good starting point to study the general 
phase diagram for the dynamics at finite temperature and finite force. On the contrary, the metric constraint is always 
a safe choice, since by definition it cannot produce an unphysical roughness strictly larger than one. 

2. Interaction between the segments of size l*[F,T) in the effective directed trap model 

For the reasons explained in details above, we thus consider that the microscopic interactions between monomers 
are given by the metric constraint. The question is then : in the effective directed trap model described in section 
MI D 3| what is the corresponding renormalized interactions between two consecutive segments of typical length 
l*(F,T) representing local quasi-equilibrated metastable states ? In the absence of more refined arguments, we feel 
that the metric constraint is actually also an appropriate renormalized interaction, since it prevents the appearance 
of an unphysical roughness £ > 1 and it allows all possible physically acceptable values £ < 1 for the roughness. 
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